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"Th I We study the r-process nucleosynthesis in neutrino-driven winds of gravitational core 
Ph| collapse SNell. Appropriate physical conditions are found for successful r-process nucle- 
Q ■ osynthesis, which meet with several features of heavy elements discovered recently in 
metal-deficient halo stars. We find also several difficulties which are not explained in the 
d '. present wind models. We discuss quests for new insights in nuclear physics, astrophysics, 
and astronomy. 

X 

H ' 1. Introduction 

! 

Stars with various masses provide a variety of production sites for intermediate-to-heavy 
mass elements. Very massive stars > lOM© culminate their evolution by supernova (SN) 
explosions which are also presumed to be most viable candidate for the still unknown 
astrophysical site of r-process nucleosynthesis. Even in the nucleosynthesis of heavy ele- 
ments, initial entropy and density at the surface of proto-neutron stars are so high that 
nuclear statistical equilibrium (NSE) favors production of abundant light nuclei. In such 
explosive circumstances of the so called hot-bubble scenario, not only heavy neutron-rich 
nuclei but also light unstable nuclei play a significant role. 

Studies of the origin of the r-process elements are critical in cosmology, too. There 
still remains potentially a serious problem that the cosmic age of the expanding Universe 
derived from cosmological parameters may be shorter than the age of the oldest globular 
clusters. Since both age estimates are subject to the uncertain cosmological distance 
scale, an independent method has long been needed. Thorium (Th) and Uranium (U), 
which are typical r-process elements and have half-lives of 14 Gyr (Th) and 4.5Gyr (^■^'^U), 
have recently been detected along with other elements in very metal-deficient stars [|l],§ . 
If we model the r-process nucleosynthesis in these first-generation stars, these long lived 
radioactive nuclei can be used as good cosmochronometers completely independent of the 
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uncertain cosmological distance scale. 

We discuss in this article that supernova explosions of very massive stars could be a 
viable site for r-process nucleosynthesis. In sect. 2, we explain how to model the neutrino- 
driven winds. The details of explosive nucelsynthesis in the hot-bubble is discussed in sect. 
3. Although our theoretical study has confirmed a successful r-process in SNell, there are 
several findings which are still difficult to explain. We discuss several quests for nuclear 
physics and astrophysics in sect. 4, and astronomy in sect. 5. 

2. Neutrino-Driven Winds in Type-II Supernovae 

Recent measurements using high- dispersion spectrographs with large Telescopes or the 
Hubble Space Telescope have made it possible to detect minute amounts of heavy elements 
in faint metal- deficient ([Fe/H] < -2) stars [|I|]. The discovery of r-process elements in these 
stars has shown that the relative abundance pattern for the mass region 120 < A is sur- 
prisingly similar to the solar system r-process abundance independent of the metallicity of 
the star. Here metallicity is defined by [Fe/H] = log[N(Fe)/N(H)] - log[N(Fe)/N(H)]0. It 
obeys the approximate relation t/10^°yr ~ lO'^^^'^l. The observed similarity strongly sug- 
gests that the r-process occurs in a single environment which is independent of progenitor 
metallicity. Massive stars with lOM© < M have a short life, ~ 10'' yr, and eventually end 
up as violent supernova explosions, ejecting material into the interstellar medium early 
on in the history of the Galaxy. However, the iron shell in SNe is excluded from being 
the r-process site because of the observed metallicity independence. 

Hot neutron stars just born in the gravitational core collapse SNell release most of 
their energy as neutrinos during the Kelvin-Helmholtz cooling phase. An intense fiux 
of neutrinos heat the material near the neutron star surface and drive matter outfiow 
(neutrino-driven winds). The entropy in these winds is so high that the NSE favors a 
plasma which consists of mainly free nucleons and alpha particles rather than composite 
nuclei like iron. The equilibrium lepton fraction. Ye, is determined by a delicate balance 
between v^ + n ^ p + e~ and z/e + p ^ + which overcomes the difference of chemical 
potential between n and p, to reach ~ 0.45. R-process nucleosynthesis occurs because 
there are plenty of free neutrons at high temperature. This is possible only if seed elements 
are produced in the correct neutron-to-seed ratio before and during the r-process. 

Although Woosley et al. demonstrated a profound possibility that the r-process could 
occur in these winds, several difficulties were subsequently identified. First, independent 
non-relativistic numerical supernova models 0] have difficulty producing the required 
entropy in the bubble, S/k ~ 400. Relativistic effects may not be enough to increase the 
entropy dramatically Second, even should the entropy be high enough, the effects 

of neutrino absorption Ve + n ^ p + e" and Ve + A{Z^ N) A{Z + 1, — 1) + e~ may 
decrease the neutron fraction during the nucleosynthesis process. As a result, a deficiency 
of free neutrons could prohibit the r-process §]. Third, if neutrinos are massive and have 
appropriate mixing parameters, energetic and are converted into due to fiavor 
mixing. This activates the + n ^ p + e~ process and results in a deficiency of free 
neutrons. 

In order to resolve these difficulties, we have studied neutrino-driven winds in a 
Schwarzschild geometry under the reasonable assumption of spherical steady-state fiow. 
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Figure 1. Nuclear reaction flow patterns in Z-N plane, and abundances in the insets at the 
time; (a) t ~ 20 ms during the a-process (upper) , and (b) t ~ 567 ms during the r-process 
(lower). Time zero refers to the time of core bounce. Highlight abundance scales to logY, 
where Y is the number fraction of each nucleus. White dots show the nuclei included in 
the reaction network code J^,]^. The neutrino-driven wind model used is for L^, = 10^^ 
ergs/s and M = 2Mq [Q. 



The parameters in the wind models are the mass of neutron star, M, and the neutrino 
luminosity, L^. The entropy per baryon, S/k, in the asymptotic regime and the expansion 
dynamic time scale, Tayn, which is defined as the duration time of the a-process when the 
temprature drops from T ^ 0.5 MeV to 0.5/e MeV, are calculated from the solution of the 
hydrodynamic equations. Then, we carried out r-process nucleosynthesis calculations in 
our wind model. We found that the general relativistic effects make Tdyn much shorter, 
although the entropy increases by about 40 % from the Newtonian value of S/k ~ 90. By 
simulating many supernova explosions, we have found some interesting conditions which 
lead to successful r-process nucleosynthesis, as to be discussed in the following sections. 
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Figure 2. R-process abundance (solid line) as a function of atomic mass number 
A compared with the solar system r-process abundance (filled circles) from Kappeler, 
T0[| . The neutrino-driven wind model used is for L^y = 10^^ ergs/s and 



Beer, & Wisshak 

M = 2Mq. The solar system r-process abundance is shown in arbitrary unit. 



3. Nucleosynthesis 

3.1. R-process and the key reaction 

Previous r-process calculations PJTl[] had shortcoming that the seed abundance dis- 



tribution was first calculated by using smaller network for light-to-intermediate mass 
elements, and then the result was connected further to another r-process network in a 
different set of the computing run. For this reasaon it was less transparent to interpret 
the whole nucleosynthesis process. This inconvenience happened because it was numeri- 
cally too heavy to run both a-process and r-process in a single network code for too huge 
number of reaction couplings among ~ 3000 isotopes. Our nucleosynthesis calculation is 
completely free from this complexity because we exploited a fully implicit single network 
code which is applied to a sequence of the whole processes of NSE - a-process - r-process. 
Shown by white dots in Fig. 1 are the isotopes included in our reaction network [|^,|[. 

Let us remind the readers that there were at least three difficulties in the previous 
theoretical studies of the r-process. The first difficulty is that an ideal, high entropy in 
the bubble S/k ~ 400 was hard to be achieved in other simulations f^-^]. 

The key to resolve this difficulty is found with the short dynamic time scale Tdyn ~ 
10 ms in our models of the neutrino-driven winds. As the initial nuclear composition of 
the relativistic plasma consists of neutrons and protons, the a-burning begins when the 
plasma temperature cools below T ~ 0.5 MeV. The ^He(aa,7)^^C reaction is too slow 
at this temperature, and an alternative nuclear reaction path ^He(Q;n, 7)^Be(a, n)^^C 
triggers explosive a-burning to produce seed elements with A ~ 100 (Fig. la). Therefore, 
the time scale for nuclear reactions is regulated by the ^He(Q;?T,, 7)^Be reaction. It is given 
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by tn = {plY^YnX{aan — Be)) . If the neutrino-driven winds fulfill the condition 
'Tdyn < ^N, then fewer seed nuclei are produced during the a-process with plenty of free 
neutrons left over when the r-process begins at T ~ 0.2 MeV. The high neutron-to-seed 
ratio, n/s ~ 100, leads to appreciable production of r-process elements, even for low 
entropy S/k ~ 130, producing both the 2nd {A ~ 130) and 3rd {A ~ 195) abundance 
peaks and the hill of rare-earth elements {A ~ 165) (Fig. lb and Fig. 2). 
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Figure 3. The same as in Fig. 2, but for the neutrino-driven wind model of = 5 x 10^^ 
ergs/s. Solid line represents the result by using the Woosley & Hoffman rate |jl2|] of the 
'^He(an, 7)^Be reaction, and long-dashed line shows the result for this rate multiplied by 
factor 2, as suggested by the recent experiment of Utsunomiya et al. |jl3 



The three body nuclear reaction cross section for ^He(a?T., 7)^Be is one of the poorly 
determined nuclear data which may alter the r-process nucleosynthesis yields. The inverse 
process has recently been studied experimentally by Utsunomiya et al. ||13[, and photo- 
disintegration cross section of ^Be has been measured with better precision than those of 
the previous experiments. Applying the principle of the detailed balance to this process, 
one can estimate the cross section for ^He(a?T,, 7)^Be. They found that the thermonuclear 
reaction rate is almost twice as big as that of Woosley and Hoffman |12| but in resonable 



agreement with the recent compilation of Angulo et al. |Tj]. However, there still remain 
several questions on the consistency between their result and electron-scattering experi- 
ments, on the contribution from the narrow resonance J'^ = 5/2~ (2.429 MeV), etc. It 
is also a theoretical challenge to understand the reaction mechanism and the resonance 
structure because two different channels, ^Be + n and ^He + a, contribute to this process. 

Therefore, we show two calculated results in Fig. 3: The solid line displays the result 
obtained by using the Woosley and Hoffman cross section |jl2|, assuming a *Be + n 
structure for ^Be. We also calculated the r-process by multiplying this cross section by 
factor of 2 (long-dashed line). This makes a drastic change in the r-process yields in 
the 3rd {A ~ 195) abundance peak. More theoretical and experimental studies of the 
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'^He(an, 7)^Be reaction are highly desired. 

3.2. Neutrino-nucleus interactions 

Neutrino interactions with nucleons and nuclei take the key to resolve the second dif- 
ficulty which was pointed out in sect. 1. The difficulty is that the effects of neutrino 
absorptions + n ^ p + e~ and z/g + A{Z,N) A{Z + l,iV — 1) + e~ during the 
a-process may induce a deficiency of free neutrons and break down the r-process condi- 
tions P]. These two types of neutrino interactions control most sensitively the electron 
fraction and the neutron fraction in a neutron-rich environment. In order to resolve this 
difficulty, we have updated the electron-type neutrino capture rates for all nuclei and 
electron-type anti-neutrino capture rate for free protons [|T^,|T^ . 



The new r-process calculation proves to be almost invariant. One can understand this 
robustness of the succesful r-process in the following way: The specific collision time for 
neutrino-nucleus interactions is given by 

where Lj^si is the individual neutrino or antineutrino luminosity in units of 10^^ ergs/s, 
ej =< Ef > / < Ei > in MeV (z = z/g, Ue, etc.), and (a^) is the cross section 
averaged over the neutrino energy spectrum. At the a-burning site of r ^ 100 km for 
L^,5i ^ 10, e^^ = 12 MeV, and (a^) ^ lO-^^cm^, r^^(r=100 km) turns out to be ^ 240 
ms. This collision time is larger than the expansion dynamic time scale; Tdyn ~ 10 ms ^ 
ri,^(r=100 km) ^ 240 ms. Because there is not enough time for z/g's to interact with n's 
in such rapidly expanding neutrino-driven wind, the neutron fraction is insensitive to the 
neutrino absorptions. (Note, however, that the opposite condition Tdyn > could apply 
to slowly expanding winds as adopted in the previous studies []3|,|^,p!T|| . In this case the 
r-process is strongly affected by the neutrino interactions [^. ) 

One might wonder if our dynamic time scale ~ 10 ms is too short for the wind to 
be heated by neutrinos. Careful comparison between proper expansion time and specific 
collision time for the neutrino heating is needed in order to answer this question. Otsuki 
et al. have found that the supernova neutrinos transfer their kinetic energy to the wind 
most effectively just above the neutron star surface at 10km <r< 20km. Therefore, one 
should refer here to the duration time for the wind to reach the a-burning site, Theat, 
rather than to Tdyn. One can estimate this heating time scale 

rt dr . 
Thcat = / — , (2) 
Jn u 

where u is the ffuid velocity of the wind. By setting the radius of the neutron star surface 
to Ti = 10 km and rj = 100 km, we get Theat ~ 30 ms. The collision time is given by 
Eq. (|T|) by setting L^,5i ^ 10, = (e^^ + epJ/2 = (12 + 22)/2 = 17 MeV, r ^10 km, and 
{cr^) ^ 10~''^cm^. Let us compare Theat and to one another: 

^ 3.4ms < Theat ~ 30ms. (3) 

We can thus conclude that there is enough time for the expanding wind to be heated by 
neutrinos even with short dynamic time scales for the a-process, Tdyn ~ 10 ms. 



7 



3.3. Neutrino oscillation 

The third difficulty discussed in sect. 1 is that the energetic massive z/^ and Vr may 
change their flavors to emerge as z/g due to the MSW effect. Although the luminosity 
of each type of neutrino Lj (z = z/g, z/g, v^, z/^, Vr, i^r) is similar, average energies are 
different from one another; e^^ = 12 MeV, ejj^ = 22 MeV, and = = 34 MeV for 
the other flavors Because of this difference, if the flavor mixing happened, newly 
converted were more energetic so that the Ug + n — > p + e~ reaction would decrease the 
neutron fraction. 

We have recently shown that the neutrino flavor oscillation destroys the r-process con- 
dition only if the mixing parameters satisfy 0.3 eV^ < Am^ [|l7l. Although recent experi- 
ments of the atmospheric neutrinos and the missing solar neutrinos have indicated much 
smaller Am^, the LSND experiment suggests remarkably larger Am^ which overlaps with 
our interesting parameter region. We should wait for more experiments. 

3.4. Extended reaction network 

The r-process is thought to proceed after the pile up of seed nuclei produced in the 
a-process at higher temperatures Tg 5 ~ 2.5. Since charged-particle reactions, which 
reassemble nucleons into a-particles and a-particles into heavier nuclei (i.e. a-process), 
are faster than the neutron-capture flow which is regulated by beta-decays, the light-mass 
neutron-rich nuclei were presumed to be unimportant. 

However, we have recently found [[18|,|19[ that even light neutron-rich nuclei progressively 
play the signiflcant roles in the production of seed nuclei. The nuclear reaction network 
used in the previous studies included only limited number of light unstable nuclei, 

3H, ^Be, S'^B, i^'i^C, i^N, 1^0, i8'20p^ 23,24^^^ therefore need to extend the 

network code so that it covers all radioactive nuclei to the neutron-drip line. We take 
the rates of charged particle reactions from those used in the Big-Bang nucleosynthesis 
calculations [p0|-23[| and the NACRE compilation |li . 



Let us briefly discuss the results of the r-process calculations, using the extended re- 
action network |]19[. At early epoch of the wind expansion, t < a few dozens ms, both 
temperature and density are so high that the charged particles interact with one another 
to proceed nucleosynthesis around the /^-stability line in the light-mass region A < 20. 
There are plenty of protons and a-particles as well as neutrons at this epoch, and the 
main reaction flow is triggered by ^He(an, 7)^Be 



^He(an, 'y)^Be{a, ny^C{n, 7)^^C(a, ny^O{n, 7)^^0(a, n)20Ne or ^^0{a, 7)20Ne..(4) 

^He(an, j^Beia, n)i2C(n, ^Y^Cin, -fY^C{n, -fY'^C... (5) 

However, at relatively later epoch when the temperature drops below T = 0.5/e MeV 
even after the a-rich freeze out, a new reaction path [|TB],|rU| starts operating: 

3H(a, 7)^Li(n, 7)^Li(a, nY^B{n, jY^B{n, -fY^B{n, ^Y^B{n, -fY^B{e-uY^C... (6) 

This new reaction flow path also takes appreciable flux of baryon number and continuously 
produces seed nuclei. The classical r-process like flow, (n,7) followed by beta decay, has 
already started from light nuclei. This is a very different result from the previous picture 
that the r-process starts from only intermediate-mass seed nuclei A 100. 
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Figure 4. Extended nuclear reaction network for explosive nucleosynthesis [|T^. For 
a-process conditions at earlier times the main reaction flows are triggered by the 
^He(Q;n, 7)^Be reaction and two flow paths (4) and (5) follow. For r-process conditions 
at later times, the third path (6), which is triggered by the ^H(a, 7)^Li(n, 7)®Li(a, n)^^B 
chain reactions, becomes equally important. See the text for more details. 



Since we do not have much information of (2n, 7) reactions, we assumed in this paper 
that the (n,7) and (7,n) reactions are in thermal equilibrium for the flrst neutron and then 
the second neutron is subsequently captured. Although inclusion of the (2n,7) reactions 
did not change drastically the flnal result, these rates are presumed to be lower limit. More 
elaborated theoretical calculation of the (2n,7) reaction has suggested that the di-neutron 
correlation may increase these cross sections. In the present calculations the yields of 



even the most neutron- rich isotopes were found to be abundant [|T^. Photodisintegration 
reactions of ^He and ^He and their electromagnetic structure also have been studied 
experimentally [0,^]. Extensive measurements of the nuclear properties of neutron- rich 
nuclei near the drip line, including ^'^'^^B and ^^'^^F, are yet to be carried out. 

There are several branching points between (n,7) and (a,n) reactions. They are at ^^C, 
^■^O, ^^Mg, etc. Experimetal studies to measure these reaction cross sections are highly 
desirable. 
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4. Quest for Nuclear Physics and Astrophysics 

We discussed the r-process nucleosynthesis in a single set of (M, L^) of the neutrino- 
driven wind in sect. 2. In realistic supernova events, however, the neutrino luminosity 
varies slowly as a function of time during Kelvin-Helmholtz cooling phase ~ 10 sec, as 
observed in SN1987A. We therefore expand the parameter space in this section in order 
to include a number of (M, L,y)-grids in the reasonable range 1.2Mq < M < 2.0Mq and 
10^°ergs/s < < 5 X 10^^ ergs/s. Calculated results of rdyn and S are displayed in Fig. 
5 0]. Shown also are two zones for which the r-process nucleosynthesis might occur so 
that the second abundance peak {A ^ 130) and the third abundance peak {A ^ 195) 
emerge from a theoretical calculation |^. In this figure, we observe that the neutrino- 
driven winds from massive and compact proto-neutron stars may produce successfully the 
r-process elements. 

In order to compare the nucleosynthesis result with the observed r-process abundance, 
the calculated yields are integrated over the mass-weighted time history, according to the 
time evolution of the neutrino luminosity 

Lj, = L^oexp(-t/r^), (7) 

where we set L^q = 5 x 10^^ ergs s~^ and t^, = 1.0 s so that the total integrated neutrino 
energy becomes 3 x 10^^ ergs to fit the observed result in SN1987A. Assuming that the 
trajectories can be described at each time by steady flow corresponding to the grids 
shown in Fig. 5, the mass- weighted r-process yields are given |0 by 

ej 

Mej = j M{t)dt, (9) 

where Yi is the number fraction of i-th nuclide, M{t) is the mass-loss rate, and Mej is the 
total mass ejected in the neutrino-driven winds from the neutron star. 

The calculated result is shown in Fig. 6. Total mass ejection is Mej = 1.5 x 1O~^M0, and 
total mass of the r-process elements Mr-process = 1.0 x IO'^^Mq. These values are compared 
with those of the benchmark study by Woosley et al. 0, which gives Mej = 3 x 10~^Mq 
and Mr_process = 1-0 X 1O~^M0. Their calculation makes larger total mass ejection because 
their highest entropy is three times bigger and the expansion time scale is about one order 
of magnitude larger than those in the present wind model. Unfortunately, their result 
turns out to be incomplete p[ for the lack of neutrino-nucleus interactions during the 
nucleosynthesis, as discussed in the previous section. We included these interactions fully 
in the present calculation. Having the same Mr-process, we can explain the absolute mass 
of the observed r-process abundance in the solar system, taking account of the frequency 
of finding SNell in the Milky Way over 10^'^ yr. 

Our result in Fig. 6 reproduces fairly well the observed second abundance peak (A ^ 
130), the hill of rare-earth elements (A ^ 165), and the third peak (A ^ 195). However, 
there are several defects, too. The biggest defect is that the abundance around A ^ 100 ~ 
120 is overproduced by, at most, factor 50. Although this factor is much smaller than 
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Figure 5. Relation between entropy per baryon, S/k, and dynamic time scale, Tdyn, for 
various combinations of the neutron star mass, 1.2Mq < M < 2.0Mq, and the neutrino 
luminosity, 10^° < < 10^^ ergs/s, taken from Otsuki et al. Solid and broken 
lines connect the same mass and luminosity. Two zones indicated by shadows satisfy the 
approximate conditions, for = 0.4, on which the successful r-process occurs ||2^ to 
make the second abundance peak around A = 130 (lower) and the third abundance peak 
around A = 195 (upper). 



the overproduction by two or more orders of magnitudes in the previous calculations 
[Tl[] , this is a common defect in the hot-bubble scenario. This overproduction is related 
apparently to the hydrodynamics of the expanding neutrino-driven winds. These nuclei 
are produced in the a-rich freeze out. One of the hints to resolve this difficulty may be the 
bypass effects due to the neutrino interactions with very neutron-rich nuclei |^ . Further 
studies are needed to suppress this unwelcome peak. 

The second defect is a shift of the third abundance peak around A ^ 195 by a couple of 
mass units. This is also a common feature found in the previous r-process calculations 
|^J|,|lll]. These elements are the beta-decay products of extremely neutron-rich unstable 
nuclei on the neutron magic N = 126. Peak position depends on the timing of freezeout 
of the r-process and subsequent beta-delayed neutron emission. Therefore, a particular 
combination of environmental evolution of neutron-number density, Nn, and temperature, 
Tg, as well as the expansion dynamic time scale, Tdyn, might match the freeze out so that 
it results in the right position of the abundance peak. 



The third defect is the deficiency of abundance right above or below the peak elements, 
i.e. at A f» 140 ~ 150 and ~ 180 ~ 190. These deficiencies seem related to yet unseen 
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effects of deformation or strucure clianges of unstable nuclei surrounding the neutron 
magic numbers N = 82 and 126. Further extensive theoretical studies and observational 
challenge to determine the masses, lives, and beta Q-values of these nuclei are highly 
desired. 

The fourth defect is the underproduction of actinoid elements, Th-U-Pu (A = 230 
~ 240), by more than one order of magnitude. The observed high abundance level of 
these nuclei (see Fig. 6) might suggest an existence of a new magic number around N 
= 150 ~ 160. Xenon ^^^Xeys and platinum ^^^Ptny are the typical r-process elements 
on the second and third abundance peaks, which are the decay products from extremely 
neutron-rich unstable nuclei with neutron magic numbers N = 82 and 126, respectively. 
From these observations we estimate that the waiting point nucleus is located by shifting 
AN fii 7 or 9 units from the peak element. Applying the same shift AN 7 or 9 to ^^^Thi42 
and ^^^Ui46, we could assume a new magic number around N = 150 ~ 160 which may 
lead to the fourth abundance peak at A = 230 ~ 240. Actually, in the very light nuclear 
systems, a new magic number N = 16 was found in careful experimental studies of 
the neutron separation energies and interaction cross sections of extremely neutron-rich 
nuclei. Since these possibilities were not taken into account in the present and previous 
calculations, the deficiency of actinoids might be improved by modernizing nuclear mass 
formula including such effects. 

Another possibility is to make actinoid elements in neutron star mergers or the mergers 
of the neutron star and black hole binaries which have extremely small lepton fraction. 
Ye < 0.2 However, these processes do not virtually produce any intermediate-mass 
nuclei including iron, which contradicts with the fact that the observed iron abundance 
is at least proportional to the r-process elements and actinoid elements over the entire 
history of Galactic evolution ~ 10^° yr. Different possibility has recently been proposed 
in Inner mass shells of the prompt supernova explosions have very small lepton 

fraction, < 0.16, and thus some actinoid elements are produced here. 

Since we discuss only material ejected from the proto-neutron star behind the shock, 
it does not make any serious problem to see the underproduction in mass region A < 
90. Most of these intermediate-mass nuclei are ejected from the exploded outer shells in 
super novae. 



5. Quest for Astronomy 

The solar system abundance is thought to consist of many different nucleosynthesis 
products in various astrophysical sites. For this reason it is not straightforward to compare 
the "pure" process to be studied theoretically with those observed experimentally in the 
solar system. However, recent detections of extremely enhanced r-process abundances 
in several metal-deficient stars have indicated that they describe surprisingly very similar 
pattern to the solar system r-process abundance which we call a "universality". We here 
proved theoretically that the neutrino- driven winds from massive and compact proto- 
neutron stars can explain this "universality". 

Let us refer to massive neutron stars. Although a large core mass 1.7M0 < M is some- 
times claimed from observational viewpoints, several established EOSs for neutron star 
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Figure 6. Final mass- weighted r-process abundance (solid line) as a function of atomic 
mass number A compared with the solar system r-process abundance (filled circles) [p!0|] . 
The neutrino-driven wind model used is for L^, = 10^^ ergs/s and M = 2Mq. The solar 
system r-process abundance is shown in arbitrary unit. 



matter are known to stabilize these massive cores. Wanajo et al. have recently pro- 
posed a formation mechanism of massive neutron star in the delayed explosion mechanism. 
Since the r-process nucleosynthesis is very sensitive to the core mass of proto-neutron 
stars it is desirable to probe the mass dependence of their yields observationally. 

Incidentally, a large dispersion in heavy element abundances of halo stars has recently 



been observed p2| , |33| . Ishimaru and Wanajo have shown in their galactic chemi- 
cal evolution model that if r-process nucleosynthesis occurs in either massive supernovae 
> 30Mq or small mass supernovae 8 — IOMq, the observed large dispersion can be well un- 
derstood. In addition, SN1994W and SN1997D are presumed to be due to 25 Mq-AOMq 



massive progenitors because of a very low Ni abundance in the ejecta [p3,3(;]. These 



massive supernovae are known to have massive iron cores > I.SMq and leave massive 



remnants 36 



It is more critical for the r-process nucleosynthesis whether the remnant is a neutron 
star or a black hole. The recent discovery ||3^ of a strong overabundance of a-elements in 
the companion in the massive black hole binary, GRO J1655-40, suggests clearly that the 
explosion was a gravitational type core-collapse, i.e. SNell. The atmosphere of the star 
orbiting around black hole or neutron star is likely to be polluted and contaminated by the 
nucleosynthesis products in supernova event associated to the formation of these compact 
objects. Therefore, if a future spectroscopic analysis of this companion star could detect 
the r-process elements and if the mass of the black hole were determined observationally, 
it would show clear evidence for the fact that the r-process could occur in SNell which 
leave remnant black holes. 

It has recently been suggested theoretically that, even should the core mass be M ~ 
IAMq which is a typical cold neutron star mass, taking proper boundary conditions of 
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the wind at the shock r ~ 10,000 km can improve crucial conditions for the r-process 
nucleosynthesis. Extensive studies are underway |^ . 



The r-process elements in the first generation stars might also make a strong impact on 
cosmology. The actinoids like U-Th-Pu are produced only in the r-process, and for their 
long lives, which are comparable to the cosmic expansion age, these elements have ever 
been used as cosmochronometers. Cosmological models infer that the cosmic expansion 
time scale is 10^° y < % < 1-5 x 10^° y, but there are several other observations which 
suggest controversially longer Galactic age than Major source of the uncertainty arises 
from the uncertain cosmological distance scale. In order to construct a reliable nuclear 
cosmochronometer, which has a great advantage being free from the distance scale, one 
has to establish that the atmosphere of extremely metal-deficient stars has for the most 
part originated from the gas of a single supernova event which produced enough amount of 
r-process elements in the early Galaxy. Extensive theoretical studies of the r-process 
have suggested that the "universality" may hold only for elements having Z = 56 ~ 80, 
while it seems slightly broken for the actinoid elements. They have suggested also that the 
Th/U chronometer is almost free from this breaking universality. These theoretical sug- 
gestions have recently been confirmed by the detections of Th along with other r-process 
elements in several new metal-poor halo giants which were observed by using the High 
Dispersion Spectrograph equipped with SUBARU Telescope, SUBARU/HDS Fur- 
ther collaboration among nuclear physics, astrophysics, astronomy, and even cosmology 
is highly desirable. 

This work was supported in part by Japan Society for Promotion of Science, and by the 
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Ministry of Education, Science, Sports and Culture of Japan, and DoE Nuclear Theory 
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